Orbital degeneracy, Hund's coupling, and band ferromagnetism: 
effective quantum parameter, suppression of quantum corrections, 

and enhanced stability 

Bhaskar Kamble and Avinash Sing 
Department of Physics, Indian Institute of Technology Kanpur - 208016 

An effective quantum parameter is obtained for the band ferromagnet in terms 
of orbital degeneracy and Hund's coupling. This quantum parameter determines, 
in analogy with \/N for the generalized Hubbard model and 1/5 for quantum spin 
systems, the strength of quantum corrections to spin stiffness and spin-wave energies. 
Quantum corrections are obtained by incorporating correlation effects in the form 
of self-energy and vertex corrections within a spin-rotationally-symmetric approach 
in which the Goldstone mode is explicitly preserved order by order. It is shown that 
even a relatively small Hund's coupling is rather efficient in strongly suppressing 
quantum corrections, especially for large A/", resulting in strongly enhanced stability 
of the ferromagnetic state. This mechanism for the enhancement of ferromagnetism 
by Hund's coupling implicitly involves a subtle interplay of lattice, dimensionality, 
band dispersion, spectral distribution, and band filling effects. 
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I. INTRODUCTION 

Experimental studies of magnetic and electronic excitations in various band ferromag- 
netic systems continue to be of strong current interest, as evidenced by intensive neutron 
scattering studies of spin-wave excitations throughout the Brillouin zone in ferromagnetic 
manganites highlighting magnon damping and anomalous zone-boundary softening,- angle 
resolved photoemission spectroscopy (ARPES) studies of iron to investigate many-body in- 
teraction between quasiparticles at the Fermi level,-"^ and spin polarized electron energy loss 
spectroscopy (SPEELS) studies of surface spin waves in ultrathin Fe films showing strong 
spin-wave softening due to reduction of exchange interaction.- Driven by recent advances in 
the resolution of experimental probes, these studies provide valuable insight into details of 
the microscopic mechanism and characteristics of band ferromagnetism. 

Realistic multi-band calculations of spin-wave dispersion using an itinerant-electron 
model, in bulk bcc Fe e.g., have so far been carried out only in the random phase approx- 
imation (RPA) owing to the complexity of the band structure.-i^ Recently, a tight-binding 
model involving 9 orbitals (4s, 4p, and 3d) per Fe atom has been used to calculate spin-wave 
dispersion in the RPA,- and electron self-energy corrections in the ferromagnetic phase of 
iron were studied in light of recent ARPES experiments on Fe. 

On the other hand, band ferromagnetism being an intrinsically strong-coupling phe- 
nomenon, spin-wave excitations in a single-band ferromagnet are strongly renormalized by 
correlation effects, as studied recently by incorporating self-energy and vertex corrections 
within a systematic inverse-degeneracy (l/A/") expansion scheme in which the spin-rotation 
symmetry of the Hamiltonian and hence the Goldstone mode^ are explicitly preserved order 
by order beyond the RPA.^ For the single-band Hubbard model, the correlation-induced 
minority-spin spectral-weight transfer was shown to result in strong spin-wave energy renor- 
malization (quantum correction), and the interplay of lattice, band dispersion, and band 
filling effects was studied on the competition between the delocalization and exchange en- 
ergy contributions to the spin stiffness, which fundamentally determines the stability of the 
ferromagnetic state in an itinerant ferromagnet.— 

So how are quantum corrections generally affected by orbital multiplicity and Hund's 
coupling? This question is of fundamental importance in view of the multi-band nature 
of transition-metal ferromagnets, but has not been addressed so far in the literature. In 
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the special case when the Hund's couphng (inter-orbital interaction) is identical to the 
intra-orbital interaction, as obtained in the generalized A/'-orbital Hubbard model,- the 
quantum corrections are simply suppressed by the inverse- degeneracy factor l/M . However, 
for arbitrary Hund's coupling, the role of orbital degeneracy on quantum corrections to 
spin-wave excitations has not been investigated so far. 

In this paper, we will extend the above Goldstone-mode-preserving approach for the 
study of correlation-induced quantum corrections to a multi-band ferromagnet with arbi- 
trary Hund's coupling. We will show that orbital multiplicity and Hund's coupling strongly 
suppress the quantum corrections and spin-wave energy renormalization in a band ferro- 
magnet. We will further show the existence of an effective quantum parameter which, in 
analogy with 1/ S for quantum spin systems and l/M for the generalized A/'-orbital Hubbard 
model, plays the role of h in effectively determining the magnitude of quantum corrections 
in a multi-band ferromagnet. 

A variety of methods have been employed to investigate the role of orbital degeneracy and 
Hund's coupling on the stability of metallic ferromagnetism, as briefly reviewed below. The 
magnetic phase diagram was calculated by finite-temperature quantum Monte Carlo simu- 
lations within the dynamical mean- field theory (DMFT), and Hund's coupling was shown to 
effectively stabilize ferromagnetism in a broad range of electron fillings even for a symmetric 
DOS.^^ It was pointed out that this stabilization is different from the mechanism based on 
asymmetric DOS, which leads to ferromagnetism in the single-band Hubbard model. 

Finite-temperature magnetism of iron and nickel was investigated using the LDA + 
DMFT approach which combines DMFT with realistic electronic structure methods, and 
many body features of the one-electron spectra and the observed magnetic moments were 
described. The Coulomb interaction energy values used were U = 2.3(3.0) eV for Fe (Ni), 
and J = 0.9 eV for both Fe and Ni, obtained from constrained LDA calculations. 

The role of lattice structure and Hund's coupling was investigated for various three- 
dimensional lattice structures within the DMFT using an improved quantum Monte Carlo 
algorithm that preserves the spin-SU(2) symmetry. It was shown that the earlier Ising-type 
DMFT calculations^- overestimate the tendencies toward ferromagnetic ordering and the 
Curie temperature. Both the lattice structure and orbital degeneracy were found to be es- 
sential for the ferromagnetism in the parameter region representing a transition metal.— 
Other numerical techniques such as exact diagonalization,— »^ density matrix renormaliza- 
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tion group,— slave boson,— and the Gutzwiller variational scheme^^ have also been em- 
ployed. 

The special case of two orbitals per site and quarter filling yields an insulating ferro- 
magnetic state with staggered orbital ordering, as suggested by the equivalence to an "an- 
tiferromagnetic" state in the pseudo-spin space of the two orbitals, and first proposed as 
a mechanism for stabilization of ferromagnetsim by Roth.^^ However, the estimated Curie 
temperature was found to be too low for transition-metal ferromagnets by a factor of 10. 
The insulating ferromagnetic state at quarter filling has also been investigated at strong 
coupling,—!^ and using the exact diagonalization method.— "^"^ 



II. TWO-ORBITAL HUBBARD MODEL WITH HUND'S COUPLING 

We consider a degenerate two-orbital Hubbard model 

H = — tiji.C'laa^j^a + '^L/3'^icr/3 + H-C) — U ^^(SjQ.SjQ, -|- Sj/j.Sj/j) — 2J ^^(SjQ,.Sj^) , (1) 

{ij),a i i 

where a and /? refer to the two degenerate orbitals at each lattice site i and Sj^ = il)\^{cr /2)ipi^j, 
are the local spin operators for the two orbitals /i = a, /5 in terms of the fermionic operators 
— ('^It/i ^^"^ Pauli matrices cr. The hopping terms Uj = t for nearest neighbours 
and t' for next-nearest neighbours. The above model includes an intra-orbital Hubbard 
interaction U, and an inter-orbital Hund's coupling J. As our objective is to investigate the 
role of Hund's coupling on quantum corrections in the orbitally-degenerate ferromagnetic 
state, we have not retained the inter-orbital density-density interaction term VoniaUi/^, also 
conventionally included in the orbital Hubbard model. This density interaction term is 
important in the context of orbital ordering in manganites,— especially near quarter filling,— 
as expected from its structural similarity with an antiparallel-spin interaction in the pseudo- 
spin {a, (3) space, which will favor antiferromagnetic orbital ordering. 

The orbital Hubbard model above has been written in an explicitly spin-rotationally- 
symmetric form. This continuous symmetry implies the existence of Goldstone modes in 
the spontaneously-broken-symmetry state, ^ which will play a central role in our study of 
correlation effects in the ferromagnetic state in which self-energy and vertex corrections are 
systematically incorporated such that the Goldstone mode is explicitly preserved order by 
order. 
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FIG. 1: Exact diagrammatic representation of the transverse spin propagator in terms of the 
irreducible particle-hole propagator. 

III. TRANSVERSE SPIN FLUCTUATIONS 



We assume a ferromagnetic ground state with polarization in the z direction, and examine 
transverse spin fluctuations representing both collective spin-wave and single-particle Stoner 
excitations. We consider the time-ordered transverse spin fluctuation propagator in this 
broken-symmetry state: 



(2) 



where the orbital indices /i,!/ — a, (3, and the fermion spin-lowering and spin-raising opera- 

The transverse spin fluctuation propagators can be expressed exactly in terms of the irre- 
ducible particle-hole propagators 0, as shown diagrammatically in Fig. 1. The corresponding 
coupled equations are: 



(3) 



where the interaction term C/^/j,/ = U for /i' — u' and C/^/j,/ = J for /i' ^ u' is shown as the 
wavy line in Fig. 1, and summation over repeated indices is implied. 

As external probes such as magnetic field or the neutron magnetic moment couple equally 
to electron magnetic moments in the two degenerate orbitals, it is physically relevant to 
consider the sum 



X ^(q,^^) = Xaa + Xa/3 (q>^) 



(4) 



of the intra-orbital and inter-orbital contributions. Indeed, it is particularly convenient to 
solve the coupled equations for the total transverse propagator, and we obtain: 

0(q,^) 



X (q,'^) 



!-([/ + J)0(q,u;) 



(5) 
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where 

0(q,u;) = 0«„(q,u;) + 0„^(q,u;) (6) 

represents the total irreducible particle-hole propagator, for which a systematic expansion 
will be discussed below. 

IV. SYSTEMATIC EXPANSION FOR THE IRREDUCIBLE PARTICLE-HOLE 

PROPAGATOR (p 

In analogy with the inverse-degeneracy (l/A/") expansion for the generalized AT-orbital 
Hubbard model,- we consider a systematic expansion: 

for the irreducible propagator 0(q, u) in orders of fluctuations. The zeroth-order first term 
(p^^^ is simply the bare particle- hole propagator, whereas the higher-order terms 0^^^ 
etc. represent quantum corrections involving self-energy and vertex corrections, as discussed 
below. 

In the inverse- degeneracy expansion scheme,- the diagrams were systematized in terms of 
the expansion parameter l/A/", with the n^'^-order term 0'-"'' involving n powers of l/A/". Thus 
the expansion parameter l/A/" played, in analogy with 1/S for quantum spin systems, the 
role of h. In the following we will show that the dimensionless factor {U"^ + J"^) / {U + (and 
a similar factor for the Af orbital case) plays the role of the expansion parameter for arbitrary 
Hund's coupling. This allows for a continuous interpolation between the single-orbital limit 
with no Hund's coupling and the generalized A/'-orbital limit with J = U. 

Random phase approximation: Retaining only the zeroth-order term 0^°^ yields the ran- 
dom phase approximation (RPA), amounting to a "classical" description of noninteracting 
spin-fluctuation modes. As the hopping term is diagonal in orbital indices, the zeroth-order 
term involves only the intra-orbital contribution: 

(q, ^) = AO (q, ^) = — F ' (8) 

where the Hartree-Fock level band energies = Ck — crA involve the exchange splitting 
2A = {U + J)m (9) 
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between the two spin bands. The superscripts +(— ) refer to particle (hole) states above 
(below) the Fermi energy ep. Here the magnetization m = 2(Sj^) is assumed to be identical 
for both orbitals fi = a, (3 in the orbitally degenerate ferromagnetic state. For the saturated 
ferromagnet, the magnetization m is equal to the particle density n for each orbital. 

At the RPA level, the two (intra- and inter-orbital) components of the transverse spin 
propagator are easily obtained by solving the coupled equations (3), and we obtain: 



[Xc 



[q,u;jjR,PA = - 



2 [l-f/+xo(q,^) l-f/-Xo(q,^) 



[Xafl (q.^)]RPA 



1 



xoiq,^j 



xoiq,^j 



[x/3/^(q.^)]RPA 
(q>^)]RPA 



(10) 

(11) 



.py^,-nr..^ 2 Ll-f/+Xo(q,^) l-f/-Xo(q,^ 
where the two effective interaction terms above are = U ± J. As seen, the propaga- 
tors involve linear combinations of two modes, which can directly be identified as in-phase 
and out-of-phase combinations with respect to the orbital degrees of freedom. These two 
modes represent gapless (acoustic) and gapped (optical) branches, as shown below. We also 
introduce kernels for these two propagators which will be used later in the expressions for 
quantum corrections: 
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(12) 



(13) 



U + J corresponds to the usual 
Goldstone mode (acoustic branch). This is expected as the exchange splitting 2A = m{U+J) 
in the xo energy denominator involves the same effective interaction f/"*". On the other hand, 
the out-of-phase mode with effective interaction U~ = U — J yields gapped excitations 
(optical branch). A typical spectral function plot is shown in Fig[3l with the inset showing 
dispersion of the acoustic and optical branches, and onset of Stoner excitations. 

Quantum Corrections: In the following we consider first-order quantum corrections to 
the irreducible particle-hole propagator (7). We consider the relatively simpler case of a 
saturated ferromagnet in which the minority-spin particle-hole processes are absent as the 
minority-spin band is pushed above the Fermi energy due to Coulomb repulsion. In this case 
the effective antiparallel-spin interactions at lowest order reduce to the bare interactions U 
and J, and the effective parallel-spin interactions reduce to a single term involving the 




FIG. 2: The first-order quantum corrections to the irreducible particle-hole propagator. 

majority-spin particle-hole bubble. Generally, these effective interactions involve a series of 
bubble diagrams, with even and odd number of bubbles, respectively. 

Diagrammatic contributions to the first-order quantum correction 0*^^^ are shown in Figure 
[21 The interaction lines are either U or J, depending on the orbitals of the connecting fermion 
lines. The external {fi = a,f3 on the right) orbital degrees of freedom are also summed over 
to include both the intra- and inter-orbital contributions as in Eq. (6). The physical meaning 
of the four diagrams has been discussed earlier.- Diagrams (a) and (d) represent corrections 
to the irreducible propagator due to self-energy corrections to the spin-| particle arising 
from spin and charge fluctuations, respectively. The shaded part in diagram (a) represents 
the propagator XJ^p^(Q,^7). Diagrams (b) and (c) represent vertex corrections, where the 
shaded part represents the kernel Fj^p^(Q,r2) introduced in Eqs. (12,13). In diagram (b), 
the opposite-spin particle-particle interaction suppresses the spin-| particle — spin-| hole 
correlation, yielding a negative correction to 0. Diagram (b) therefore represents suppression 
of the magnetic response due to particle-particle correlations. All four diagrams involve a 
spin-charge coupling, as indicated by the spin-| particle-hole bubble, present explicitly in 
diagrams (c) and (d) and imphcitly in (a) and (b). 
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Integrating out the fermion frequency-momentum degrees of freedom, the first-order 
quantum corrections to the irreducible particle-hole propagator are obtained as: 

{iu' + J')x-JiQ,n) + 2UJx-JiQ,n)} 

2 
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-k"-q+Q "^k" 

We note that in the J ^ limit of decoupled orbitals, we recover the single-band Hubbard 
model results.- 

As collective spin- wave excitations are represented by poles in (5), spin- rotation symmetry 
requires that (j) = 1/{U + J) for g, = 0, corresponding to the Goldstone mode. Since the 
zeroth-order term (j)^^^ already yields exactly !/([/ + J), the sum of the remaining terms 
must exactly vanish in order to preserve the Goldstone mode. For this cancellation to hold 
for arbitrary J and U, each higher-order term 0*^"^ in the expansion (7) must individually 
vanish, implying that spin-rotation symmetry is preserved order-by-order. This cancellation 
is demonstrated below for the first-order quantum correction obtained above. 
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Towards this end, we first note that the boson term (quantity in braces) in Eq. (15) for 
(j)^''^ can be expressed as: 

t/r-+(Q, n) + JT-jiQ, n) = {iu' + j')x-^ + 2UJx-J}/xo (is) 

which is of identical form as the boson terms in Eq. (14) and (16) for 0*^"^ and (f)^^\ With 
ej^^q — ej^~ = 2A for g = 0, we obtain from Eqs. (14-17): 

0(1) (q = 0, cu) = (f)^''^ + (j)^'^ + (j)^'^ + (f)^''^ 

-2 {{U' + j')x-: + 2UJx-;} 
+ {iU' + J')iXa^-Xo) + 2UJx-;} 

+ {iU' + J')xo}] , (19) 

which yields identically vanishing contribution for each spin-fluctuation mode Q, thus pre- 
serving the Goldstone mode. We note that this mode-by-mode cancellation is quite indepen- 
dent of the spectral-weight distribution of the spin-fluctuation spectrum between collective 
spin-wave and particle-hole Stoner excitations. Furthermore, the cancellation holds for all 
uj, indicating no spin-wave amplitude renormalization, as expected for the saturated ferro- 
magnet in which there are no quantum corrections to magnetization. 

Evaluation of the Q integral in Eqs. (14-17) has been discussed earlier.— Using a spectral 
representation it is convenient to carry out the Q integral numerically so as to include all 
three (acoustic, optical, and Stoner) contributions from the magnon [x ^,r~"'") and the 
particle-hole terms. A typical spectral function plot of the RPA-level magnon propagator 
(Q) '^)]rpa shows (Fig. [3]) the low-energy (acoustic), intermediate-energy (optical), and 
high-energy (Stoner) contributions, with the inset showing the dispersion of acoustic and 
optical branches and onset energy of the Stoner branch. 

V. HUND'S COUPLING AND SUPPRESSION OF QUANTUM CORRECTIONS 

We will now show that Hund's coupling results in a strong suppression of quantum cor- 
rections, and this is the central message of this paper. The suppression is approximately 
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FIG. 3: All three contributions — the low-energy (acoustic), intermediate-energy (optical), and 
high-energy (Stoner) — as seen in a typical spectral function plot of the magnon propagator 

[XaaiQ^ ^)]rpa for Q = (§, §), are included in evaluating the Q integral in Eqs. (14-17). Inset 

1 

shows the dispersion of the acoustic and optical branches and the onset energy for the Stoner 
branch, for the sc lattice with t' = 0.25, U = 15, J = 3, and band filling n = 0.5. 

by an overall factor of (f/^ + J'^)/{U + J)^ for the two orbital case, and by a factor of 
{U^ + {Af - l)j2)/(t/ + (A/" - l)Jf for the J\f orbital case. The magnitude of the overall 
factor is asymptotically exact in the two limits J and J U . In these two limits the 
factor approaches 1 and 1/2, respectively, for the two orbital case, whereas it approaches 
1 and l/A/", respectively, for the Af orbital case. Indeed, for our degenerate two-orbital 
Hubbard model (1), this overall factor plays the same role as l/A/" in the generalized M- 
orbital Hubbard model,- and generalizes this earlier result to the case of arbitrary Hund's 
coupling. In the generalized A/'-orbital Hubbard model, the form —U(Yli^Si^).{Yli^Siy) of 
the interaction term implies identical inter- and intra-orbital interactions {J = U). 

To see this overall suppression factor, we consider Eqs. (14-17) for the quantum correc- 
tions (f). In Eq. (14), assuming similar contributions from the acoustic and optical modes, 
which is a good approximation in the J U limit when the two modes become nearly 
degenerate, the boson term 2U JXap yields negligible contribution in view of the opposite 
contributions (11) of the acoustic and optical modes. Hence, only the contribution from 
the first boson term {U"^ + J'^)Xaa survives, leaving an overall factor {U"^ + J^) on carry- 
ing out the Q, VL integration. Comparing with the corresponding factor {U + J)^ for the 
single orbital case with identical interaction {U + J), yields an overall relative factor of 
(t/^ + J'^)/{U + jy. In view of Eq. (18), the quantum correction 0^'') also involves the same 
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boson term and hence yields the same overall factor. Similarly, the quantum corrections 0^*^^ 
and ^^'^^ together yield the same boson term and hence the same overall factor. 

The magnitude of the overall factor is also exact in the J ^ U limit, as shown below. 
The boson term in Eq. (14) can be identically written, in view of Eqs. (10,11), as: 

iU' + J')X-J + 2UJx~J = l[(U + Jfx-Jo + iU- J)^Xopt] (20) 

in terms of the acoustic and optical branches. Therefore, in the J ^ U limit, the contribution 
of the optical mode vanishes, leaving an overall factor of 1/2, as also resulting from the 
expression {U'^ + J'^)/{U + Jf. 

Band ferromagnetism being a strong-coupling phenonomenon, quantum corrections to 
spin stiffness and magnon energy for the single-orbital Hubbard model with U ^ W have 
been shown to yield strong renormalizations (reduction relative to the RPA values) in two 
and three dimensions.— The strong suppression of quantum corrections shown above high- 
lights the critical role of Hund's coupling in stabilizing ferromagnetism in realistic systems 
like transition metals such as Fe, Ni etc. with multiple "id orbitals. We propose that in such 
systems, Hund's coupling favours ferromagnetism by strongly suppressing the quantum cor- 
rections. 



VI. QUANTUM CORRECTIONS TO SPIN STIFFNESS 

We now consider the net quantum correction <^ for small q in order to obtain the renor- 
malized spin stiffness, which provides a quantitative measure of the stability of the ferro- 
magnetic state with respect to long- wavelength fluctuations. In Eqs. (14-17), writing the 
antiparallel-spin particle-hole energy denominator as 

eic-q - = 2A[1 + (ek-q - ek)/2A] (21) 
and expanding in powers of the small band-energy difference 

6 = -(ek_q - ek) = q.Vek - ^{q.Vfe^ , (22) 

we find that besides the zeroth-order cancellation for q = 0, the first-order terms in S 
also exactly cancel. This exact cancellation implies that there is no quantum correction 
to the delocalization contribution (V^ek) in the spin stiffness constant; only the exchange 
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FIG. 4: Quantum correction to spin stiffness for the simple cubic lattice, shown as a function of 
band filling for different J with fixed U + J = 1.5W = I8t. 

contribution in the spin stiffness is renormalized by the surviving second-order terms in S, 
and we obtain for the first-order quantum correction to stiffness: 
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(23) 



where d is the lattice dimensionality. Here we have introduced effective interactions in the 
transverse channel: 



u^M,^) = {iu' + j')x-^iQ,n) + 2UJx-;iQ,m 
u:M,n) = u:M,n)~{u' + j')x\Q,n). 



(24) 



As Vek is odd in momentum, the second and third terms in (23) give vanishingly small 
contributions due to partial cancellation. 

The behaviour of quantum correction to spin stiffness with band filling is shown in Fig. 
4 for different values of the Hund's coupling J. Here we have kept U + J fixed so that 
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FIG. 5: Rapid suppression of quantum correction to spin stiffness with Hund's coupling J, sliown 
for the sc lattice at two different band fillings. Inset shows comparison of the normalized quantum 
correction with the approximate form (1 + (J/C/)^)/(l + J/U)"^, which asymptotically approaches 
the calculated result in the two limits J /U and J/U 1. 

the exchange band splitting and the classical spin stiffness remain unchanged. The quan- 
tum correction to spin stiffness decreases sharply with J and reduces to exactly half the 
magnitude when J = U, as shown in Fig. 5 for two different band fillings. Inset shows a 
comparison of the normalized quantum correction calculated from (23) with the approximate 
form (1 + (J/f/)^)/(l + J/U)'^, which asymptotically approaches the calculated result in the 
two limits J/U ^ and J/U 1, as discussed in section V. These results clearly show 
that even a small Hund's coupling is rather efficient in strongly suppressing the quantum 
correction. 

Fig. 6 shows that the renormalized spin stiffness D = D^'^^ — D^^'^ rapidly changes 
sign from negative to positive with increasing J, highlighting the effective stabilization of 
ferromagnetism by Hund's coupling, shown here for the simple cubic lattice with t' = 0.25 
and fixed U + J= 1.5W. 

As discussed in section V, generalizing to the A/'-orbital case, the quantum correction to 
spin stiffness should be approximately suppressed by the factor {U"^ + {Af— 1) J^)/(?7+ (A/" — 
1) J)^. With increasing Hund's coupling, this factor rapidly approaches l/A/", particularly for 
large A/"; the generalized A/'-orbital Hubbard model therefore provides a good approximation 
for transition-metal ferromagnets with A/" = 5 3d orbitals. 
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FIG. 6: The effective stabilization of the ferromagnetic state by Hund's coupHng is shown by the 
rapid change with J from negative to positive spin stiffness, shown as a function of band fihing for 
the simple cubic lattice with t' = 0.25 and fixed U + J = LbW . 
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FIG. 7: Renormalized spin stiffness for different number of orbitals A/", showing the \/J\f suppression 
of quantum corrections with orbital degeneracy, evaluated for the bcc lattice with bandwidth 
W = 16t = 3.2eV, Coulomb interaction energy U = W = 3.2eV, and lattice parameter a = 2.87A 
for Fe. The measured value for Fe is 280 meV.A^. 

We have examined the role of this l/A/" suppression of quantum corrections on the spin 
stiffness. Fig. 7 shows the renormahzed spin stiffness D = D^^^ — jfD^^^ for different number 
of orbitals A/", evaluated for a bcc lattice with t'/t = 0.5, bandwidth W = 16t = 3.2eV, 
Coulomb interaction energy U = W = 3.2eV, and the lattice parameter a = 2.87A for Fe. 
In a recent band-structure calculation,- the interaction energy considered is U = 2.13eV (so 
that the magnetic moment evaluated per Fe atom is equal to 2.12/iB) and the bandwidth from 
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FIG. 8: Hund's coupling results in a strongly momentum-dependent enhancement of magnon 
energies, as seen from the magnon dispersion along symmetry directions in the Brillouin zone for 
the sc lattice at two different band fillings with t' = 0.25 and fixed U + J = 18t = 1.5W. The 
rapid crossover from negative- to positive-energy long-wavelength modes with J shows the strong 
stabilization of ferromagnetism by Hund's coupling. 



the calculated DOS plot is seen to be about 4eV. Our calculated values for the renormalized 
spin stiffness for A/" = 5 are close to the measured value 280meVA^ for Fe. The spin stiffness 
is seen to involve a quantum reduction of about 25% near optimal filling. 



VII. RENORMALIZED MAGNON DISPERSION 



Turning now to the magnon dispersion over the entire Brillouin zone, the renormalized 
magnon energy ciJq for finite q is obtained from the pole condition 1 — (f/+ J)Re 0(q, — c^^q) = 
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FIG. 9: Rapid suppression of quantum correction and stabilization of ferromagnetism due to Hund's 
coupling, as seen in the magnon dispersion along symmetry directions for the square lattice at two 
different band fillings with t' = 0.5 and fixed U + J = 12t = 1.5W. The dispersion along the T-X 
direction shows pronounced anomalous softening near the zone boundary. 

in Eq. (5), where 0(q, cj) = </)*^°^(q, tu) + (f)^^\q,uj), and the four contributions to the first- 
order quantum correction (f)^^\q,uj) are given in Eqs. (14-17). While the bare particle-hole 
propagator (p^^^q^, to) remains real in the relevant cu range, the quantum correction 0*^^^(q, to) 
is complex for any finite u < due to the coupling with charge fluctuations, resulting in finite 
zero-temperature magnon damping. Both collective and Stoner excitations are included in 
evaluating the Q integral, as discussed below (19). 

The effect of Hund's coupling on the renormalized magnon dispersion is shown in Figs. 8 
and 9, which provide comparisons of the bare (tUq) and renormalized (cjq) magnon dispersions 
along symmetry directions for the simple cubic and square lattices, respectively. Again U+J 
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was kept fixed so that the bare magnon energy remains unchanged. While Hund's couphng 
increases the magnon energy in the entire Brillouin zone, the effect is particularly dramatic 
for long wavelength modes, which rapidly crossover from negative-energy to positive-energy 
modes. 

Even when the bare magnon dispersion exhibits nearly Heisenberg form, with energies at 
X,M,R approximately in the ratio 1:2:3 as in Fig. 8(b), the renormalized magnon dispersion 
shows strong anomalous softening at X relative to R. This indicates that magnon renor- 
malization due to spin-charge coupling results in the "generation" of additional exchange 
couplings J2, J3, J4 etc. within an equivalent localized-spin model with the same magnon 
dispersion. 

The renormalized spin stiffness and magnon energies obtained above essentially deter- 
mine the finite-temperature spin dynamics and therefore the Curie temperature Tg. Finite- 
temperature spin dynamics in a band ferromagnet and reduction of magnetization due to 
thermal excitation of spin waves has been recently discussed in terms of the spectral weight 
transfer across the Fermi energy.— Due to the spin-flip scattering of electrons accompanying 
the thermal magnon excitations, a portion of the majority-spin spectral weight is transferred 
to the minority-spin band above the Fermi energy, while an equal amount of minority-spin 
spectral weight is transferred to the majority-spin band below the Fermi energy. From the 
self-energy correction in Eq. (14), and using the resolution of the transverse spin propaga- 
tors (Eqs. 10,11) into acoustic and optical modes, the Curie temperature is approximately 
given by: 

1 1 



q ^ q qj k Vk-q 



(25) 



within a renormalized spin-fluctuation theory,— where cUq and cj* refer to the magnon en- 
ergies for the acoustic and optical modes. Now, finite Hund's coupling results in i) gapped 
optical modes with higher energies (cj* > cj*) than the acoustic modes, ii) reduced weightage 
of the optical modes due to the {U — J)^ factor in above equation, and iii) enhancement 
of the acoustic magnon- mode energies (even at fixed U + J). All three factors suppress 
the rhs of above equation, directly resulting in enhanced Curie temperature. While this 
enhancement of Tc is in general agreement with DMFT studies where only local excitations 
are incorporated, Eq. (25) highlights the sensitivity to long-wavelength modes as well in our 
Goldstone-mode-preserving approach. Negative-energy long-wavelength modes as in Figs. 8 
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and 9 would result in vanishing even if bulk of the (short-wavelength) modes have posi- 
tive energy. Long wavelength modes play a particularly important role in low-dimensional 
systems. The divergence of the rhs of Eq. (25) in one and two dimensions yields vanishing 
Tc in accordance with the Mermin- Wagner theorem.— 

VIII. CONCLUSIONS 

The role of orbital degeneracy and Hund's coupling on quantum corrections to spin-wave 
excitations in a band ferromagnet was investigated. A spin-rotationally-sjTiimetric approach 
was employed in which self-energy and vertex corrections are incorporated systematically 
so that the Goldstone mode is explicitly preserved order by order. The present study of 
quantum corrections for arbitrary Hund's coupling allows for a continuous interpolation 
between the orbitally independent case ( J = 0) equivalent to the single-band Hubbard model 
and the orbitally symmetric case of identical inter- and intra-orbital Coulomb interactions 
{J = U) equivalent to the generalized Hubbard model with Af degenerate orbitals per site, 
for which the first-order quantum corrections are suppressed by the factor I/JV. We find that 
even a relatively small Hund's coupling is rather efficient in strongly suppressing the quantum 
corrections, especially for large A/", resulting in a strong enhancement of ferromagnet ism. 

This mechanism for the enhancement of ferromagnet ism due to strong suppression of 
quantum corrections by Hund's coupling implicitly involves an interplay of several band 
and lattice characteristics. Competition between the delocalization (V^ek) and exchange 
((Ve^)^) /2A contributions to spin stiffness results in a subtle interplay of lattice, dimension- 
ality, band dispersion, spectral distribution. Coulomb interaction, and band filling effects, 
as investigated for several two and three-dimensional lattices.— 

For a two-orbital Hubbard model, the first-order quantum correction to spin stiffness 
was obtained diagrammatically and evaluated as a function of the Hund's coupling strength 
J. We showed that the effective quantum correction factor decreases rapidly from 1 at 
J = to 1/2 at J = ?7, and its behaviour with J closely follows the approximate form 
{U^ + J^)/ {U + jy obtained from our diagrammatic analysis. In the intermediate-coupling 
regime, with interaction strength U comparable to the bandwidth, the renormalized spin 
stiffness was evaluated in the whole range of band fillings, and was found to rapidly crossover 
from negative to positive values with increasing Hund's coupling, highlighting the strong 
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role of orbital degeneracy and Hund's coupling on the stability of the ferromagnetic state 
with respect to long wavelength fluctuations. We also obtained, for both the square and 
simple cubic lattices, the renormalized spin-wave energy dispersion cjq for momenta along 
symmetry directions in the Brillouin zone. While Hund's coupling results in an enhancement 
of the magnon energy in the entire Brillouin zone, the effect is particularly dramatic for long 
wavelength modes, which rapidly crossover from negative to positive energy. 

Generalizing to the A/" orbital case, an effective quantum parameter (t/^ + (A/'— 1) J^)/(f/+ 
(A/"— was obtained which, in analogy with 1/S for quantum spin systems and l/A/" for 
the generalized Hubbard model, plays the role of h for quantum corrections in a band ferro- 
magnet. For large Af, this quantum parameter decreases rapidly with the Hund's coupling 
J and saturates to l/A/" as J — ?7. 

This strong suppression of quantum corrections due to orbital degeneracy and Hund's 
coupling is quite significant for the 3d transition-metal ferromagnets Fe, Co, Ni, where 

= 5. With J/U = 1/4, as considered in the recent RPA calculations for iron,- we obtain 
a value of 5/16 ~ 0.3 for the quantum parameter. While the smallness of this quantum 
correction parameter accounts for why RPA calculations of the spin stiffness for Fe with 
realistic band structure yield values in close agreement with the measured value of 280 
meV.A^, it also highlights the significant magnitude of the correlation-induced quantum 
corrections involved in the measured spin stiffness values for transition-metal ferromagnets. 
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